Ionization heating in rare-gas clusters under intense XUV laser pulses 
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The interaction of intense extreme ultraviolet (XUV) laser pulses 

(A = 32 nm, / = 10^^~^ W/cm^) with small rare-gas clusters (Ari47) is studied by quasi- classical 
molecular dynamics simulations. Our analysis supports a very general picture of the charging and 
heating dynamics in finite samples under short- wavelength radiation that is of relevance for several 
applications of free-electron lasers. First, up to a certain photon flux, ionization proceeds as a 
series of direct photoemission events producing a jellium-like cluster potential and a characteristic 
plateau in the photoelectron spectrum as observed in [Bostedt et aL, Phys. Rev. Lett. 100, 
013401 (2008)]. Second, beyond the onset of photoelectron trapping, nanoplasma formation leads 
to evaporative electron emission with a characteristic thermal tail in the electron spectrum. A 
detailed analysis of this transition is presented. Third, in contrast to the behavior in the infrared 
or low vacuum ultraviolet range, the nanoplasma energy capture proceeds via ionization heating, 
i.e., inner photoionization of localized electrons, whereas collisional heating of conduction electrons 
is negligible up to high laser intensities. A direct consequence of the latter is a surprising evolution 
of the mean energy of emitted electrons as function of laser intensity. 

PACS numbers: 52.50.Jm, 36.40.Wa,41.60.Cr 



I. INTRODUCTION 

The rapid advances in producing ultrashort and in- 
tense short-wavelength radiation keep stimulating to- 
day's laser-matter research. Whereas sub-fs vacuum ul- 
traviolet (VUV) pulses generated from ultrashort near 
infrared (NIR) laser fields have become a key tool in at- 
tosecond physics |1|, coherent and highly intense pulses 
from free electron lasers (PEL) are currently pushing the 
frontiers of intense laser-matter science from the VUV 
up to the x-ray domain ^^. Promising examples of 
novel applications that are expected to become feasible 
with upcoming x-ray FEL's are diffr active or holographic 
imaging of finite samples with unprecedented spatial and 
temporal resolution [^|-g| . A key fundamental issue being 
of relevance for many applications of short-wavelength 
PEL light is how the nature and dynamics of intense 
laser-matter interactions develop as function of wave- 
length, e.g., from the infrared towards the x-ray domain. 

Promising prototype systems for corresponding ex- 
ploratory studies on finite targets are atomic clusters, 
as great experience in their experimental and theoretical 
investigation is available [|l0|-|l2|. Their efficient absorp- 
tion in the well-studied NIR range (A ^ 800 nm) results 
in energetic reaction products, such as, highly charged 
ions (l3[ 0, fast electrons g [§, x-rays @ |l|, or 
even neutrons [^. In this domain the response is gov- 
erned by a strongly confined and short-lived nanoplasma, 
whose formation is typically triggered by field ioniza- 
tion. Particularly strong absorption occurs for resonant 
collective excitation of the nanoplasma electrons. In- 
tense laser-cluster experiments at shorter wavelengths. 
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i.e., in the VUV (100 nm), became feasible by turn of 
the millennium at the DESY free electron laser (FLASH) 
and have shown unexpectedly strong heating and highly 
charged ion production [|o|, though collective plasma 
heating effects become negligible due to low quiver am- 
plitudes of electrons. In particular, strong polarization 
fields and resonance enhancement as typical for the NIR 
disappear [^ . A further fundamental difference of the 
VUV domain is that vertical photoionization becomes 
the trigger ionization process for the nanoplasma for- 
mation. As electronic screening and local fields from 
neighboring ions in the cluster lower the atomic ioniza- 
tion barriers, cluster constituents can be charged much 
stronger than isolated atoms with the same photon en- 
ergy via photoionization |2^, |23| . Taking this effect into 
account, the observed cluster heating could be explained 
by inverse Bremsstrahlung (IBS) in the resulting dense 
nanoplasma [^. Additional aspects that have been put 
forward to explain strong absorption in the VUV in- 
clude the effects of realistic ion potentials for describing 
IBS [Q and enhanced heating through many-body colli- 
sions in a transient strongly coupled nanoplasma [^ . A 
general conclusion b eing in line with the observed ther- 
mal electron spectra [^ is that laser heating of quasifree 
electrons via IBS is a leading absorption process in the 
A ~ 100 nm regime. Simulations on ArN support that 
IBS remains substantial even with A = 60 nm p3| . 

Evidence for a transition to a new response regime 
have been found in experiments in the extreme ultra- 
violet range (XUV) [g^^ where photoelectron spectra 
from rare-gas clusters (ArN, (N) = 80) under 30 fs FEL 
pulses at A = 32 nm showed pronounced nonthermal, 
plateau-like signatures. Up to laser intensities of about 
lO^^W/cm^ the experimental electron spectra can be 
well explained with a multistep ionization model that ex- 
cludes nanoplasma effects and assumes sequential atomic 



photoionization in the developing cluster Coulomb poten- 
tial. An additional low energy tail emerging for higher 
intensities could only partly be explained with the mul- 
tistep model, suggesting a dynamical mixture of direct 
photoemission and thermal electron evaporation. 

The purpose of the current study is to examine the 
intensity-dependent ionization and heating processes at 
A = 32 nm {hw =38 eV) in more detail. A quasi-classical 
molecular dynamics approach including plasma effects 
like atomic ionization threshold lowering through local 
plasma fields is applied to the model system Ari47. The 
simulation results indicate that the energy capture of the 
nanoplasma at high fiuences is dominated by a new mech- 
anism we term ionization heating. Therein, the absorp- 
tion is no longer mediated by laser heating of nanoplasma 
electrons, but by inner photoionization of remaining lo- 
calized electrons. A clear marker for this heating mech- 
anism in the electron emission spectra is discussed. 

The remainder of the text is organized as follows. Me- 
thodical details of the simulations are subject of Sec. p. 
Section. |l| proceeds with a presentation of numerical 
results and discusses characteristic phases of the interac- 
tion process, the intensity dependence of electron energy 
spectra, and signatures from ionization heating. Conclu- 
sions are drawn in Sec. IV. 



II. METHODS 

As a full quantum mechanical description of clusters 
in strong laser fields is currently out of reach, quasi- 
classical approaches are widely used to model the in- 
teraction process, see, e.g., Refs. [^ Q and references 
therein. The common strategy is to describe atomic ion- 
ization processes quantum mechanically via appropriate 
rates, whereas the dynamics of the resulting ions and 
electrons is treated classically. The classical propagation 
of electrons is justified for sufficiently strong excitations 
and has been applied with great success to laser-cluster 
interactions from the NIR up to the x-ray domain. We 
thus follow this general strategy in our simulation code. 

Starting point of a simulation run is the initialization 
of an Ari47 cluster in relaxed icosahedral structure. For 
the annealing a binary Lennard- Jones potential has been 
used |2^. At this stage only neutral atoms appear ex- 
plicitly in the code. As optical field ionization is negligi- 
ble in the XUV range, only atomic photoionization is in- 
cluded as inner ionization process. We assume sequential 
one-photon absorption to be the dominant channel and 
neglect multi-photon ionization. Ionization events are 
sampled statistically in each timestep At using a Monte- 
Carlo scheme with the probability for the k-th atom/ion 
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where Wk = ^^o-'^ I{t)/JiijJ gives its actual ionization 
rate for the instantaneous laser intensity I{t) = Iq f'^{t) 
summed over the relevant atomic levels a (3s and 



TABLE I: Level occupation (occ), ionization energies (/p), 
and photoionization cross sections (a) at huo = 38 eV of 
atomic Ar as function of charge state. Electron removal from 
the 3s or the 3p level is considered (as indicated). Values 
with superscripts have been taken from *NIST, "[^, and 
^Ipd]. Remaining ionization potentials are calculated with an 
atomic all-electron dirac-LDA code ||3l|]. 



3p levels are included here) with the correspond- 
ing photoionization cross sections <j^. Here /q and 
f{t) = exp(— 21n2t^/r^) denote the peak intensity and 
the normalized temporal field envelope of a gaussian laser 
pulse with pulse duration r (FWHM), respectively. 

The photoionization cross sections are taken from free 
atomic Ar, see Tab. |, for which Ar^^ is the maximum 
charge state that can be attained via one-photon pho- 
toionization with huj = 38 eV. In a cluster, however, fur- 
ther ionization is possible because of the presence of the 
local plasma environment |2^, |3| . Considering an atomic 
ion in the charged cluster, neighboring ions and electronic 
screening by quasifree cluster electrons lead to an effec- 
tive ionization threshold / * = /p — Aenv that is lowered 
with respect to the pure atomic value by the environ- 
mental shift Aenv Here /* specifies the minimum energy 
required to lift the electron into the quasi-continuum of 
the cluster. The shift Agnv is evaluated directly from 
the plasma field in the MD simulation using the scheme 
developed in Ref. [^ . Therefore the local ionization bar- 
rier for a given ionic cell is assumed to be located at half 
distance to the next neighbor (rbar)- The shift Agnv is 
then given by the electron potential produced from all 
charges (including the residual ion) at the point of the 
barrier (rbar) minus the electron potential in the envi- 
ronmental field alone (residual ion removed) sampled at 
the position of the ion. Ionization is tested if the con- 
dition hu > /* is fulfilled. Cross sections for ionization 
steps that are energetically forbidden in free atomic Ar 
are approximated by the values for the highest possible 
atomic ionization step scaled with the remaining num- 
ber of electrons in the tested atomic level. For electron 
removal from the 3s level we assume instantaneous relax- 
ation of the residual ion to the ground state. However, 
because of a much larger cross section, 3p photoioniza- 
tion is by far the dominant process in our case. 

If an ionization event occurs, the atomic charge state 
is incremented and an electron is born at the position of 
the residual ion. The electron momentum is determined 
from energy conservation and initialized with random ori- 
entation. Active particles, i.e., all plasma electrons and 
ions, are propagated classically under the infiuence of 



the mutual Coulomb interaction and the laser field using 
a standard Verlet-algorithm. The coupled equations of 
motion read 



miVi = ^ieEia 
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where tti^, g^e, and r^ are the mass, charge, and position 
of the i-th particle and Eias = ezEQf{t) cos{ujt) is the 
electric field of a linearly polarized laser pulse in dipole 
approximation with peak amplitude Eq = y^zlo/cso. 
Here c and Eq denote the vacuum speed of light and the 
vacuum permittivity. The pairwise Coulomb interaction 
Vij is described with a pseudopotential of the form 
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with the elementary charge e, the interparticle distance 
rij , their charge states qi and qj , and a numerical smooth- 
ing parameter s. The latter regularizes the Coulomb in- 
teraction and offers a simple route to avoid classical re- 
combination of electrons below the lowest possible quan- 
tum mechanical energy level. We use a constant value 
of s and choose it to be the minimal value, such that 
the binding energy of an electron to a g-fold charged ion 
is always lower than the q-th ionization potential. For 
the particular case of argon this yields s = 1.128 A [ p3| . 
We like to note that unphysical recombination occuring 
for lower values of the smoothing parameter can lead to 
artificial heating. 

The above described MD scheme includes the full clas- 
sical dynamics of plasma electrons and ions in the clus- 
ter and resolves thermalization and IBS heating in the 
nanoplasma. To analyze their impacts on the electron 
emission we will compare the MD simulations to re- 
sults obtained from the simpler multistep model [ p7| . 
The latter neglects any nanoplasma effects and assumes 
that cluster ionization proceeds as a series of instanta- 
neous electron emission events due to direct photoemis- 
sion from the developing cluster Coulomb field. Ioniza- 
tion is treated in the same way (via Monte-Carlo sam- 
pling) as in the MD code, but ionization events are ac- 
cepted only if the single particle energy of the released 
electron is positive. The electrons are immediately re- 
moved from the simulation by neglecting any further en- 
ergy exchange (instantaneous escape). In the final spec- 
trum they occur with the single particle energy 
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at the instant of birth (electron removal from level a 
of the i-th ion). The last term in Eq. (H) describes the 
Coulomb downshift due to previously generated ions with 
charge states qj at positions Vj . For brevity, the resulting 
spectra are henceforwards termed Monte-Carlo results 
(MC) as opposed to the fully dynamical MD approach. 



III. RESULTS AND DISCUSSION 

In the following we present numerical results for the 
excitation of Ari47 by 30 fs pulses at A = 32 nm, i.e., for 
parameters similar to the experiments of Ref. [^. To 
obtain sufficient statistics, ensemble averaging over 10^ 
independent simulations is performed. Hence, all observ- 
ables reflect statistical mean values. 



A. Phases of the Ionization Dynamics 

To provide a basis for the forthcoming discussions, we 
begin with the analysis of characteristic stages of the in- 
teraction process. Fig. ^ displays the evolution of selected 
observables as predicted by MC and MD simulations for 
two representative laser intensities. For MD both the 
number of innerionized electrons nu (equal to the total 
number of photoionization events) as well as the number 
of continuum electrons ncont (positive single particle en- 
ergy, equal to outer ionization) are displayed. For MC 
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FIG. 1: (color online) Calculated ionization dynamics of Ari47 
for excitation with 30 fs laser pulses at A = 32 nm for two laser 
intensities (as indicated). The evolutions of total inner ioniza- 
tion riii and the number of continuum electrons ricont are given 
for MD and MC simulations. For the lower intensity case all 
lines lie on top of each other. The shaded area depicts the 
laser intensity envelope. Dotted horizontal lines indicate ion- 
ization levels for partial (^par) and full (^fuii) frustration (see 
text). The inset displays the evolution of the average kinetic 
energy (^kin) of trapped electrons (negative single particle 
energy) from MD runs for I — 2x lO^^W/cm^ with fixed and 
free ionic structure, respectively. 



only the latter is given and serves as a reference for clus- 
ter charging when nanoplasma effects are suppressed. 

For the lower intensity of Iq = 2 x 10^^ W/cm^, outer 
ionization evolves identical for both the MD and the MC 
approach. Electrons are excited directly to the contin- 
uum and practically no active electrons remain bound. 
This supports the claim of direct quasi-instantaneous 
electron emission without notable collision induced elec- 
tron trapping. Hence, outer ionization is nearly propor- 
tional to the integrated pulse intensity and images the 
number of photons absorbed via photoionization. The 
dynamics can be fully described within the multistep pic- 
ture. 

The ionization behavior changes substantially for the 
higher intensity of /q = 2 x 10^^ W/cm^. While MC and 
MD results still match in the rising edge of the pulse, the 
predictions begin to deviate at about t = — 15fs. From 
this point on, inner ionization starts to exceed outer ion- 
ization (MD). This reflects the onset of space charge trap- 
ping, where direct photoemission becomes partially frus- 
trated, i.e., nanoplasma generation begins. 

For a homogenous charge distribution in the cluster 
(jellium approximation) the trapping potential reads 
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where q is the total cluster charge (outer ioniza- 
tion prior to the considered ionization event) and 
R = VgN'^^^ is the cluster radius determined from the 
atomic Wigner-Seitz radius (r^^ = 2.2lA). Partial frus- 
tration of direct photoemission begins in the cluster 
center if l^eii(<7par, ^ = 0) + ficj — /p < , i.e., as soon 
as electrons are born with negative single particle en- 
ergy. Considering 3p photoionization of a neutral Ar 
atom in the cluster, partial frustration sets in for a 
cluster charge state ^par = 12, in good agreement with 
the data in Fig. |[ Subsequent ionization steps from 
regions where direct emission remains operational in- 
crease the depth of the space charge potential and 
eventually induce full frustration of direct photoemis- 
sion if l^eii(^fuii, r = R) -\- hu — Ip < 0. Note that the 
thresholds for partial and full frustration are related by 
Qiuii/Qpar = 3/2 within this electrostatic approximation. 
In our case, full frustration is reached for a cluster charge 
of ^fuii = 18, which reflects the upper limit for multistep 
ionization. 

Whereas the MC data shows smooth convergence of 
outer ionization to the full frustration limit gfuii, a tran- 
sient saturation at a slightly lower value is found near 
t = — 5 fs in the corresponding MD result. This indicates 
a small enhancement of the trapping field by electrons 
that are produced in the center of the cluster, i.e., with 
insufficient energy to escape. More importantly, a second 
long-living emission feature emerges around t ~ fs and 
outer ionization finally increases to values well beyond 
g'fuii. This additional emission persists after the laser 
pulse and can be attributed to thermal electron evap- 
oration from the nanoplasma. 



For a closer analysis of this evaporative emission we 
examine the evolution of the nanoplasma in more de- 
tail. The total number of trapped electrons is given 
by ^trapped = ^ii - ^cont, cf. Fig. |l|. Notc that the 
nanoplasma population strongly exceeds the number of 
directly emitted electrons for the higher intensity case. 
The maximum evaporation rate of thermal electrons oc- 
curs near t = 15 fs, i.e, even before the maximum number 
of trapped electrons is reached. The reason for that is the 
competition between population and temperature in the 
nanoplasma. The evolution of the average kinetic en- 
ergy of trapped electrons shows a heating phase up to 
about t = 15 fs and rapid cooling afterwards, see inset of 
Fig. |l|. The peak of this average thermal energy well co- 
incides with the maximum evaporation rate, underlining 
the strong temperature dependence of evaporative emis- 
sion. Two processes contribute to the nanoplasma cool- 
ing: (i) the evaporation itself and (ii) cluster expansion. 
The latter is driven by both Coulomb forces from clus- 
ter charging as well as thermal pressure and results in 
very efficient expansion cooling of trapped electrons. A 
comparison of runs with free vs. frozen ions shows that 
expansion cooling is the dominant process, see inset of 
Fig. 0. The rapid decay of the evaporation rate after the 
laser pulse is therefore mainly due to expansion cooling 
of nanoplasma electrons. Note that thermal evaporation 
is almost completed around t = 150 fs. 

The above analysis shows that the XUV induced elec- 
tron emission proceeds in two fundamentally different 
stages, i.e., first via multistep ionization and subse- 
quently by evaporative emission from the nanoplasma. 
For reaching the evaporative regime, the cluster ioniza- 
tion has to be sufficiently high to frustrate direct photoe- 
mission, i.e., to generate the nanoplasma. 



B. Electron Emission Spectra 

In the next step we analyze characteristic signatures 
of the above described ionization stages in the electron 
emission spectra. Fig. || displays the evolution of the final 
energy distributions of emitted electrons as function of 
laser intensity (from left to right, upper panels) for MC 
and MD simulations. An energy vs. birth time analysis 
(middle and lower panels) allows assignment of particular 
features to different electron release times. 

For intensities up to about / = 5 x lO^^W/cm^, the 
respective MC and MD spectra are almost identical, see 
Figs. ||(a)-||(c), and reflect the typical behavior of multi- 
step ionization. The multistep process produces a shoul- 
der at the low energy side of the atomic photoline, which 
broadens and develops into a plateau as intensity in- 
creases. Note that the main contribution is due to the 
emission of 3p-electrons, cf. Fig. ^(a). Starting at the 
energy of the atomic photoline E = hw — Ip (first elec- 
tron), electron energies are continuously downshifted for 
each ionization step due to the cluster space charge, as 
reflected by the signatures in the energy vs birth time 
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FIG. 2: (color online) Intensity dependent electron spectra from Ari47 (from left to right). The upper panels (a-f) compare 
final energy distributions calculated from MC (dotted) and MD (solid) simulations. Middle and lower panels display an energy 
vs. birth time analysis. Therefore the final energy distributions (color-coded) and mean energies (solid curve) of electrons are 
given as function of their release time from the mother ion. Hereby the middle row [panels (g)-(l)] displays the data from MC 
runs whereas the lower row reflects corresponding MD results. Characteristic signatures from multistep ionization and thermal 
electron evaporation are indicated in panels (q) and (r). The laser intensity envelope is schematically shown (dashed). 



evolutions, see Figs. ||(g)-||(i) and Figs. ^m)-||(o). Note 
that the latter are practically equal for MC and MD (mid- 
dle and lower panels) up to / = 5 x lO^^W/cm^. The 
emergence of the plateau in the energy spectrum can be 
understood by recalling the trapping potential in jellium 
approximation given in Eq. (||). According to this, emit- 
ted electrons show an average energy downshift with re- 
spect to the atomic photoline of 
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which scales linearly with charge state. As the average 
energy distribution for a given ionization step exhibits a 
finite width (except for the first electron), the contribu- 
tions add up to a plateau-like spectrum eventually. 

For intermediate laser intensity [Fig. ||(d)] an addi- 
tional feature appears at low electron energies on top 
of the plateau in the MC and the MD spectra. It can be 
traced back to selective photoemission resulting from par- 
tial frustration, where only electrons from outer regions 
of the cluster can be emitted directly. The resulting re- 
duced energy downshift per ionization step increases the 
signal at low energies. Full frustration of direct pho- 
toemission at higher intensities is reflected in the rapid 
saturation of MC spectra, see Figs. ||(d)-||(f). Though 
the MD spectra begin to substantially exceed the MC re- 
sults in this intensity range due to thermal evaporation, 
the characteristic energy downshift with birth-time (mul- 



tistep feature) persists in ear lyst ages for all scenarios, see 
Figs. |(j)|(l) and Figs. ip)|](r). 

The departure from pure multistep behavior for 
/ > lO^^W/cm^ is in agreement with experimental find- 
ings [^ and numerical results from kinetic Boltzmann 
simulations [Q . The roughly exponential additional con- 
tribution in the energy spectra due to electron evapora- 
tion exhibits a characteristic signature in the birth-time 
analysis. Comparison of the energy vs birth time plots 
of MD and MC runs reveals that thermal evaporation 
shows up as an additional broad feature in the MD re- 
sults (see Figs. g(p)-Fig. ^r) as soon as full frustration 
is reached. The onset of evaporation induces a turnover 
of the average electron energy vs birth-time (solid curve) 
from a decaying to an almost constant evolution. The 
latter effect is most pronounced in Figs. |^(q) and ||(r). 
Interestingly, the average final energy of thermally emit- 
ted electrons is only weakly dependent on the laser in- 
tensity. Focussing on the time after the pulse peak, the 
average energy increases by less than a factor of two when 
going from 2 x 10^^ to lO^^^W/cm^, cf. Figs. |(q)-|(r). 
Though more intense pulses produce a higher number 
of thermal electrons, their average energy changes only 
weakly. This is a central result of our analysis and indi- 
cates a characteris tic he ating process. We come back to 
this aspect in Sec. IIIC. 



Both the persistence of multistep ionization at high 
intensity (in early stages) and the emergence of evapo- 




FIG. 3: (color online) Energy correlation analysis of emitted 
electrons for the MD simulations of Fig. |2|. The plots dis- 
plays the final kinetic energy distributions (^fin) of electrons 
as function of their single particle energies at the instant of 
birth (^birth)- Characteristic features from multistep ioniza- 
tion and thermal evaporation are indicated. 



rative emission are corroborated by an analysis of the 
correlation of single particle energies of electrons at 
the instant of birth (£^birth) with their final energy af- 
ter emission (^fin). Figure contains a correspond- 
ing correlation plot for the MD runs of Fig. g. Up 
to about / = 5 X 10^^ W/cm^ the birth- and final en- 
ergies are highly correlated (almost equal), forming a 
characteristic multistep line around En^i = £^birth > 0, 
cf. Figs. ||(a)-||(c). The latter is found in all runs, but 
its increasing width at higher laser intensity indicates 
the growing effect of energy exchange between multi- 
step electrons during their escape. Emitted electrons 
with initial birth energies £^birth < begin to appear for 
/ > 5 X 10^^ W/cm^. Their uncorrelated Efi^ and £^birth 
values reflect the (random) energy exchange with other 
nanoplasma electrons required for evaporative emission. 
Besides being a marker for electron emission processes, 
the energy correlation analysis contains an image of the 
cluster potential during the interaction time with the 
laser pulse. The lower £^birth cutoff of the distribution 
is a measure for the depth of the trapping field in the 
cluster. With increasing intensity the minimum birth en- 
ergies decrease and reach values of about £^{JJJ^] 



birth 



-90 eV 
at / = lO^'^W/cm^ see Fig. |(f). The fact that the depth 
of the space charge potential becomes much larger than 
the photon energy reflects the growing contribution of 
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FIG. 4: (color online) Intensity dependent heating and ion- 
ization of Ari47 in 30 fs laser pulses at A = 32 nm. The upper 
and middle panel show the relative fraction of IBS heating to 
the total energy absorption and the average energy of emit- 
ted electrons, respectively. The lower panel displays the final 
values of total inner ionization (n\s) and the number of contin- 
uum (emitted) electrons (ricont) for MD and MC simulations. 
Dashed horizontal lines indicate ionization levels for partial 
(^par) and full (^fuii) frustration. Vertical dotted lines indicate 
laser intensities of the scenarios discussed in Figs. and pi 



temperature-driven charge separation to the cluster po- 
tential (thermal spill-out of nanoplasma electrons). 



C. Signatures of ionization heating 

We now return to a closer analysis of the interrelation 
of ionization and heating processes in the nanoplasma. 
Figure |4|(c) displays a systematic comparison of final in- 
ner (n^i ) and outer cluster ionization (p^^x) froi^ MD 
runs with outer ionization from MC simulations (n^^^) 
as function of intensity. The similar evolutions of outer 
ionization from MC and MD reflect pure multistep ion- 
ization up to about /o = 3 X 10^^ W/cm^. Note that the 
partial frustration threshold is approximately reached for 
this intensity. Beyond that intensity the nanoplasma for- 
mation begins (deviation of inner and outer ionization in 
the MD results), but multistep ionization prevails until 
full frustration sets in near /q = lO^^W/cm^. Up to 
this point the mean energy of emitted electrons (MD), 
see Fig. §(b), decreases steadily with laser intensity and 
approaches a value of 50% of the energy of the atomic 
3p photoline, reflecting roughly the center of mass of the 
multistep plateau in the energy spectra, cf. Fig. g. For 
higher intensity thermal evaporation [difference between 



outer ionization of MD and MC in Fig. ^(c)] becomes sub- 
stantial and exceeds the fraction of multistep electrons 
near Iq = 8 x lO^^W/cm^. The degree of inner ioniza- 
tion grows steadily with pulse power and reaches values 
up to {qa) = 5.5 (per atom) in the displayed intensity 
range. This results in a very dense nanoplasma and un- 
derlines the importance of ionization threshold lowering 
in the cluster. 

However, even if a dense nanoplasma is formed, the rel- 
ative impact of IBS heating to the total energy absorption 
(MD) remains negligible (< 1% for the displayed inten- 
sity range), see Fig. Q(a). This is in strong contrast to 
the VUV domain, where nanoplasma heating via IBS is 
substantial and typically even the dominating absorption 
channel [^ ^. The vanishingly small contribution of 
IBS heating in the XUV case implies that the excess en- 
ergy from inner photoionization acts as the main source 
of thermal energy in the nanoplasma. We thus suggest to 
term this energy absorption process ionization heating. A 
direct consequence of the latter is that the nanoplasma 
temperature cannot be pushed unhmitedly by just in- 
creasing the laser intensity. Because of the temperature 
dependence of electron evaporation, heating processes 
leave clear traces in the intensity dependent energy of 
emitted electrons, see Fig. ||(b). In fact, the presence 
of ionization heating leads to a very unusual evolution 
of the mean electron energy, showing an almost constant 
value from the onset of nanoplasma generation up to high 
intensities. The excitation with intense short- wavelength 
pulses thus offers a route to produce transient nanoplas- 
mas with well controlled density (tunable by the inten- 
sity) at nearly constant temperature (adjustable with the 
photon energy). 



IV. SUMMARY AND CONCLUSIONS 



intensities - at least in the rising edge of the pulse - has 
been demonstrated and estimates for the thresholds for 
partial and full frustration of direct photoemission have 
been given based on a simplified jellium picture. Ther- 
mal electron evaporation occurs as soon as a nanoplasma 
is formed. The latter requires frustration of direct pho- 
toemission and is thus delayed by a certain number of 
ionization steps. The major heating mechanism of the 
nanoplasma is a process we term ionization heating^ i.e., 
inner photoionization of localized electrons, while IBS 
absorption of plasma electrons is negligible. 

Both multistep ionization and ionization heating are 
generic mechanisms that are of relevance for a wide class 
of laser-particle interactions in the wavelength range from 
the XUV up to x-ray domain. For example, intense laser 
irradiation required for single-shot x-ray diffractive imag- 
ing of finite samples produces a deep trapping poten- 
tial due to multistep ionization, which is crucial for sec- 
ondary ionization processes. A sufficiently strong trap- 
ping field triggers field ionization cascades, which can 
lead to charge migration processes |3^ and substantial 
depletion of the number of bound electrons. A suppres- 
sion of the latter is of importance for retrieving use- 
ful diffraction pattern from the sample [p6|. Another 
direction is the controlled generation of small plasma 
targets, whose density and temperature could in prin- 
ciple be tuned over a wide range. While the intensity 
mainly controls the density (if the frustration threshold 
is reached), the photon energy can be varied to control 
the temperature produced via ionization heating. We 
thus expect that the mechanisms of multistep ionization 
and ionization heating will not only be of fundamental 
importance for laser-cluster science, but have far reach- 
ing implications for the rapidly growing field of intense 
short-wavelength laser-matter interactions in general. 



The stages and intensity dependence of charging and 
heating in rage-gas clusters under intense XUV laser 
pulses have been analyzed by comparison of molecular 
dynamics simulations with statistical Monte-Carlo re- 
sults. We find that the electron emission proceeds in 
two sequential phases: first via direct photoemission and 
then via thermal electron evaporation. The applicability 
of the concept of multistep ionization up to high laser 



AcknoAvledgments 

We like to thank Christian Gnodtke and Ulf Saalmann 
for inspiring discussions with regard to the interpreta- 
tion of the birth time analysis. Financial support by the 
Deutsche Forschungsgemeinschaft within the SFB 652/2 
and computer time provided by the HLRN Computing 
Center are gratefully acknowledged. 



[1] F. Krausz and M. Ivanov, Rev. Mod. Phys. 81, 163 

(2009). 
[2] V. Ayvazyan et a/., Eur. Phys. J. D 37, 297 (2006). 
[3] W. Ackermann et a/., Nature Photon. 1, 336 (2007). 
[4] Stanford Linac Coherent Light Source (LCLS), USA, 

http://lcls.slac.stanford.edu. 
[5] European XFEL project, http://www.xfel.eu/de/. 
[6] SPring-8 Compact SASE Source (SCSS), Japan, 

http://www.riken.jp/XFEL/eng. 
[7] R. Neutze, R. Wouts, D. van der Spoel, E. Weckert, and 



J. Hajdu, Nature 406, 752 (2000). 
[8] H. N. Chapman et a/.. Nature 448, 676 (2007). 
[9] H. N. Chapman et a/.. Nature Phys. 2, 839 (2006). 
[10] V. P. Krainov and M.B. Smirnov, Phys. Rep. 370, 237 

(2002). 
[11] U. Saalmann, C. Siedschlag, and J. M. Rost, J. Phys. B 

39, R39 (2006). 
[12] Th. Fennel, K.-H. Meiwes-Broer, J. Tiggesbaumker, P.- 

G. Reinhard, P. M. Dinh, and E. Suraud, Rev. Mod. 

Phys. 82, 1793 (2010). 



[13] M. Lezius, S. Dobosz, D. Normand, and M. Schmidt, 

Phys. Rev. Lett. 80, 261 (1998). 
[14] E. M. Snyder, S. A. Buzza, and A. W. Castleman Jr., 

Phys. Rev. Lett. 77, 3347 (1996). 
[15] E. Springate, S. A. Aseyev, S. Zamith, and M. J. J. 

Vrakking, Phys. Rev. A 68, 053201 (2003). 
[16] V. Kumarappan, M. Krishnamurthy, and D. Mathur, 

Phys. Rev. A 67, 043204 (2003). 
[17] A. McPherson, B. D. Thompson, A. B. Borisov, K. Boyer, 

and C. K. Rhodes, Nature 370, 631 (1994). 
[18] Ch. Prigent, C. Deiss, E. Lamour, J. P. Rozet, D. Vern- 

het, and J. Burgdorfer, Phys. Rev. A 78, 053201 (2008). 
[19] T. Ditmire, J. Zweiback, V. P. Yanovsky, T. E. Cowan, 

G. Hays, and K. B. Wharton, Nature 398, 489 (1999). 
[20] H. Wabnitz et a/.. Nature 420, 482 (2002). 
[21] D. Bauer, J. Phys. B 37, 3085 (2004). 
[22] Ch. Siedschlag and J. M. Rost, Phys. Rev. Lett., 93, 

043402 (2004). 
[23] L Georgescu, U. Saalmann, and J. M. Rost, Phys. Rev. 

A 76, 043203 (2007). 
[24] R. Santra and C. H. Greene, Phys. Rev. Lett. 91, 233401 

(2003). 
[25] Ch. Jungreuthmayer, L. Ramunno, J. Zanghelhni, and 

Th. Brabec, J. Phys. B. 38, 3029 (2005). 
[26] T. Laarmann, M. Rusek, H. Wabnitz, J. Schulz, A. R. B. 

de Castro, P. Giirtler, W. Laasch, and T. Moller, Phys. 

Rev. Lett. 95, 063402 (2005). 
[27] C. Bostedt, H. Thomas, M. Hoener, E. Eremina, T. Fen- 
nel, K.-H. Meiwes-Broer, H. Wabnitz, M. Kuhlmann, 



E. Plonjes, K. Tiedtke, R. Treusch, J. Feldhaus, A. R. B. 
de Castro, and T. Moller, Phys. Rev. Lett. 100, 133401 
(2008). 

[28] J. F. Ogilvie and F. Y. H. Wang, J. Mol. Struct. 273, 
277 (1992). 

[29] U. Becker and D. A. Shirley, VUV and Soft X-Ray Pho- 
toionization (Springer, 1996). 

[30] A.M. Covington et aL, XXII International Conference 
on Photonic, Electronic and Atomic Collisions, Santa Fe, 
New Mexico, U.S.A., Abstracts of Contributed Papers, 
p. 48 (2001). 

[31] A. L. Ankudinov, S. I. Zabinsky, and J.J. Rehr, Comp. 
Phys. Comm. 98, 359 (1996). 

[32] Th. Fennel, L. Ramunno, and Th. Brabec, Phys. Rev. 
Lett. 99, 233401 (2007). 

[33] The same smoothing parameter for electron-electron, ion- 
ion, and electron-ion interactions is used to avoid force 
discontinuities for surrounding particles during an ion- 
ization event. A natural but often acceptable drawback 
of this treatment is that the effect of energetic collisions 
at low impact parameters is underestimated. 

[34] B. Ziaja, T. Laarmann, H. Wabnitz, F. Wang, E. Weck- 
ert, C. Bostedt, and T. Moller, New J. Phys. 11, 10 
(2009). 

[35] Ch. Gnodtke, U. Saalmann, and J. M. Rost Phys. Rev. 
A 79, 041201(R) (2009). 

[36] T. Nakamura, Y. Fukuda, and Y. Kishimoto, Phys. Rev. 
A 80, 053202 (2009). 



